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Abstract 

In this paper, we consider a Diophantine quasi-periodic time-dependent analytic per¬ 
turbation of a convex integrable Hamiltonian system, and we prove a result of stability of 
the action variables for an exponentially long interval of time. This extends known results 
for periodic time-dependent perturbations, and partly solves a long standing conjecture of 
Chirikov and Lochak. We also obtain improved stability estimates close to resonances or 
far away from resonances, and a more general result without any Diophantine condition. 


1 Introduction and results 

1.1 Introduction 

Let n > 1 be an integer, D C M” an open bounded convex domain and T” := M"/(27rZ)". 
Consider a smooth Hamiltonian function H defined on the domain T” x H of the form 


H{e,I) = h{I)+ef{9,I), e>0, (0, /) = (0i,...,A,...,4) G T” x L>, (1) 

and its associated Hamiltonian system 

\h{t) = -dg^H{e{t),I{t)) = -edeJ{9{t)J{t)) “ “ 

For e = 0, the system is stable in the sense that the action variables I{t) of all solutions 
are constant, and all solutions are quasi-periodic. Now for e ^ 0 but sufficiently small, a 
fundamental result of Nekhoroshev states that if the system is real-analytic and the integrable 
part h satisfies a generic condition (called steepness), then the action variables I{t) of all 
solutions are almost constant for an interval of time which is exponentially long with respect 
to the inverse of the perturbation. More precisely, the following estimates hold true along all 
solutions: 
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for some positive constants Ro,TQ,c,a and b. We refer to |Nek77| and [Nek79] . The most 
important constants appearing in the estimates ([5]) are undoubtedly the constants o and b, 
which are called the stability exponents. The simplest class of steep integrable Hamiltonians 
are the convex (or quasi-convex) Hamiltonians, and these exponents depend then only on 
the number of degrees of freedom n. Nekhoroshev’s original proof yielded the following 
dependence 

, 1 

a ~ 0 ~ —;r. 

Such values for the exponents were however much worse than the values 



conjectured by Chirikov f [Chi79| i on a basis of a heuristic argument and numerical simu¬ 
lations. This issue was later solved by Lochak f |Loc92j i: more precisely, Lochak-Neishtadt 
f [LN92| i and independently Poschel ( |Pbs93 ] ) proved that ([2]) holds true with the values 


a = b = 


1 

2n 


More generally, it is proved in [BMll] that for any 0 < (5 < (2n(n — 1)) 


one can choose 


1 

“ 2(n- 1) 


S, b = 5(n — 1), 


recovering the latter result by letting 5 = (2n(n — 1)) Examples of Arnold diffusion 
m es 96], [B es97] . [Zhallj i show that in any event a < (2(n —2)) ^ so that in the convex case, 
the estimate ([2|) is, as far as the dependence on a is concerned, quite sharp. 

Now instead of an autonomous perturbation as in ([T|), one may consider a periodic time- 
dependent perturbation, that is one looks at the Hamiltonian 


H{e,I) = h{I) e>0, {e,I)eT^xB, t£T. (3) 

Setting t = tp G T and introducing a variable J E M canonically conjugated to ip, it is 
equivalent to study the autonomous Hamiltonian 

H{9,I,ip,J) = h{I) + J + ef{e,I,ip), e>0, (6,1) gT^ x B, {ip,J)GTxR. (4) 

Indeed, {9{t),I{t),ip{t),J{t)) is a solution of the system associated to (jH) if, and only if, 

{9{t),I{t)) is a solution of the system associated to ([3|). The Hamiltonian (jH) has now re -|- 1 

degrees of freedom, the integrable part h{I) + J is no longer convex but it is quasi-convex: thus 

the results of |LN92] and [P bs93] apply, and the estimates ([2|) hold true for the Hamiltonian ([3|) 

with ^ 

a = b = -r. 

2(re -I- 1) 

Now a periodic time-dependent perturbation is nothing but a special case of a quasi-periodic 
time-dependent perturbation, and we may more generally consider a Hamiltonian of the form 


H{9,1) = h{I) + £f{9, /, to), e > 0, 


{9,I)gT^xB, = G'I^ 


(5) 
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where a E is a vector which is assumed to be non-resonant, that is A: • a 7 ^ 0 for any 
non-zero k E and where • denotes the Euclidean scalar product. We will assume actually 
that 0 . satisfies a Diophantine condition: there exist 7 > 0 and r > m — 1 such that 


|A:-a|>7|A:| A; = (fci, ..., E Z”'\ {0}, |A:| := |A:i| H- \km\. (Dio-y,.r) 

As before, setting (p = ta ^ T™ and introducing a vector J = (Ji,..., Jm) E M”* canonically 
conjugated to p, the Hamiltonian ([5]) is equivalent to 


H{e,I,p,J) = h{I) + a-J + ef{ej,p), e > 0 , (0,/)eT"xH, (vp, J) e T™ x M’”. ( 6 ) 


The Hamiltonian ([ 6 ]), in the special case m = 1, reduces to Hamiltonian (jH): indeed, a E M 
and by a scaling one may assume that a = 1 , and moreover (Dio^y^,-) is obviously satished for 
7 = |q!| = 1 and r = m — 1 = 0 . 

We can now state the most general form of a conjecture of Chirikov (see [Chi79] . |CV89] 
and [CV96| 1. stated in a more precise manner by Lochak f |LMS03] l. 


Conjecture 1. The estimates m hold true for the Hamiltonian (l6|), provided it is real- 
analytic and h convex, with the exponents 


a = b 


1 

2 (n -b 1 + r)' 


Note that this conjecture is made plausible by the fact that in the periodic case, that is 
m = 1 and r = 0, it is a theorem. Yet for m > 2 and hence r > 1, it is an open question 
whether estimates ([2]) hold true for the Hamiltonian Q, for some values of a and b whatsoever. 

It is the purpose of this article to solve this problem: we will prove that the estimates ([2]) 
hold true for the Hamiltonian ([6]), provided it is real-analytic and h convex, with the exponents 
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We refer to Theorem [1] below for a more precise statement. Concerning the values of the 
exponents we obtain, let us just make two comments (a more detailed discussion is contained 
in Section First, for the periodic case (m = 1, r = 0), we also recover the known values 
of the exponents, so our result can also be considered as a “correct” generalization of the 
periodic case. Then, in the non-periodic case (m > 2, r > 1), our exponent a is substantially 
worse than the one of Conjecture [U but at the same time our exponent b is always much 
better: b is close to 1/2 (it is always strictly bigger than 1/4), and when n is hxed and m 
(or r) becomes arbitrary large, it becomes arbitrarily close to 1/2. Moreover, in the non¬ 
periodic case, b is essentially independent of n which is an interesting feature. In any case, 
it is therefore still an open question whether the values of the exponents of Conjecture [T] can 
be reached or not. 

In fact, we do obtain more general results. First, the perturbation will be allowed to 
depend also one the J variables, that is we can replace f{6,I,p) in ([ 6 ]) by f{9,I,p,J), 
provided / is bounded and real-analytic in J, when J varies in This will be the precise 
content of Theorem [TJ To understand the interest of this seemingly mild extension, one can 
compare such a result with the preservation of invariant tori (that is, the KAM theory) for 
Hamiltonians as in ([H) , Q or ([U]) . In the autonomous case of ([I]) , if the integrable Hamiltonian 
is convex then it is in particular Kolmogorov non-degenerate and the classical KAM theory 
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applies (see |Pos01| for a survey). In the periodic case ([H), the integrable Hamiltonian is no 
longer Kolmogorov non-degenerate but as it is quasi-convex, it is Arnold (or iso-energetically) 
non-degenerate: it follows that tori are preserved at a fixed energy for (j 3 |) yielding invariant 
tori for Q. In this case the perturbation may also depend on J € M without affecting 
the result. Now in the more general quasi-periodic case ([ 6 ]), the integrable part is both 
Kolmogorov and Arnold degenerate. However, using the non-degeneracy with respect to the 
I variables and the fact that the perturbation is independent of J, it is not hard to prove, 
using classical KAM techniques, that many tori with prescribed Diophantine frequencies of 
the form (w, a) E M”'"*'™' are preserved (up to our knowledge, this was first observed by 
Galavotti in |Gal94| in a restricted situation and later by Lochak in |Loc98j . |LMS03j in a 
general situation). It is crucial here to have a perturbation which is independent of J: if not, 
the method simply breaks down and it is rather easy to construct counter-examples (as in 
|Sev03j l to the preservation of (full dimensional) invariant tori. Therefore unlike the situation 
in KAM theory, our result is exactly the same when / is allowed to depend on the J variables. 

Then, as in the autonomous or periodic case, we do obtain enhanced stability close to 
resonances. In our situation, the multiplicity of any resonance is at most d, where 0 < d < n 
(with the convention that any frequency is resonant of multiplicity at least 0 ), and solutions 
who start sufficiently close to such a resonance satisfy ([ 2 ]) with the exponents 
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In the case of a resonance of maximal multiplicity d = n these exponents read 


a(n) 
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2((n + l)r + l)’ 


6(n) 


1 

2 ' 


This will be proved in Theorem [2j One recovers the improved stability exponents of the 
periodic case by setting r = 0 , and our main result by setting d = 0 . 

The complement of the neighborhoods of all resonances is the non-resonant domain (this 
domain contains, in particular, invariant tori, if any). Solutions starting in the non-resonant 
domain satisfy ([ 21 ) with the exponents 
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Moreover, the complement of this non-resonant domain is actually very small: its measure is 
of order 

K jn -b 1 )t -b 1 

2 (n + l)(l + r) 

and therefore goes to zero with e. This will be stated as Theorem [3j This actually improves 
on the corresponding statement for r = 0 , where the measure estimate of the complement is 
just of order one. This non-resonant domain comes from the proof of our main theorem, and 
ends up quite large as we need to exclude very small neighborhoods of resonances. A more 
natural definition of a non-resonant domain yield the following result: given any 0 < 7 ' < 7 
and any t' such that t' > n + m — 1 and t' > r, where 7 and r are the constants appearing 
in (Dio.y^T-)) there is a set whose complement has a measure of order 7 ', such that on this set, 
the estimates ([ 21 ) hold true with the exponents 
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This will be the content of Theorem 01 Hence on a smaller non-resonant subset, but which is 
still relatively large, we have a stronger stability result, with a' arbitrarily close to the value 
conjectured and h' much better. As a matter of fact, when r > m — 1 (for m > 2, the set 
of vectors a for which r = m — 1 has zero measure), one can choose t' = n + t and then a' 
coincides with the value conjectured. 

Finally, as usual with results in Hamiltonian perturbation theory concerning long but 
finite time scale, the Diophantine condition (Dio.-j.^T-) on the vector a G turns out to be 
unnecessary. For an arbitrary vector a G which is assume to be non-resonant, that is 


k-a^O, 

we will obtain in Theorem 0] a more general stability result which reduces to the main result 
in the case where a is Diophantine. As a matter of fact, even the assumption that a is non¬ 
resonant is unnecessary: indeed, a being fixed, if it is resonant we can write, without loss of 
generality, a = {a, 0) G x for some 1 < k < m — 1, with a G R^ non-resonant 

and the above applies immediately with a G R™ replaced by d G R^. 


1.2 Main result 

Let us now state more precisely our main result. We consider a Hamiltonian of the form 

= h{I) + a- J + (0,/)gT”xD, ((^, J) G T™ X R"* (H) 

where h is the integrable part and / the perturbation. The Hamiltonian h, defined on D := 
D X R™ by 

h{I,J):=h{I)+a-J, {I,J)eD 

will be called the extended integrable part. The functions h and / are assumed to be real- 
analytic as follows. Given two parameters tq > 0 and sq > 0, we dehne the complex domains 

Vr,D:={IeC^ I ||/-D|| <ro}, :={(/, J) G | ||(/, J) - D|| < tq} 

where 

||/-D||:= inf ||/-/'||, ||(/,J)-D|| := inf _ ||(/, J) - J')!!, 

and 


V^jn+m ^ C”+”*/(27rZ)”+”* 


max |lm(0j)| < sq, max |Im((/?j)| < sq}. 

l<i<n l<i<m 


Let us also define the associated real domains 

UroD := VroD n R”, Ur^D := Vr^D n R”+”^ = Ur^D X R"*. 

The function h is assumed to be real-analytic on Vr^D so that h is real-analytic on Vr^D, and 
its Hessian V^/i is assumed to be uniformly bounded on the complex domain Vr^D, namely 
there exists M > 0 such that 

sup \\V^h{I)\\= sup \\V^h{I,J)\\ < M (M) 

{I,J)eVroD 
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where the matrix norm is the one induced by the Euclidean norm. The gradient of h is also 
assumed to be uniformly bounded on the real domain Ur^D, that is there exists > 0 such 
that 

sup ||(v/i(/),a)|| = sup ||v/i(/)|| < n. (n) 

i&Ur^D leUr^D 

Moreover, the integrable Hamiltonian is assumed to be (strictly, uniformly) convex: there 
exists m > 0 such that for any v € M”", 

V‘^h{I)v ■ V > m\\v\\‘^. (m) 


Observe that m < M, and, without loss of generality, we may assume that m < 1. 

Finally, the function / is real-analytic on Vr^D x and moreover, given a small 

parameter e > 0, it is assumed that 

\f\ro,so<£ (e) 

where the Fourier norm \ f\ro,so of / is defined as follows: letting 


be the Fourier expansion of / with respect to {0,(p), we dehne 

|/|,„,,„:= sup _ \{k,l)\ = \k\ + \l\. 

(I,J)&VrQD 

We can now state our main theorem. 

Theorem 1. Let H be as in dSl), with h satisfying pTl) . (fHI) and (jmj), and f satisfying (j^. 
Assume also that a satisfies (Dio.y^T-); o,nd let us define 
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If e < min{eo,e*}, for any solution {I{t), J{t),9{t),ip{t)) of the system associated to H with 
initial condition (Iq, ^Oj V^o) £ D x we have 


\\I{t)-Io\\<R{e):=R.(j-^ , |t|<r(£):=T,exp(|(^)“). 

Let us just make one simple comment concerning the assumption (j^, which requires 
the perturbation f{9, ip, /, J) to be uniformly bounded in J, with J belonging to the vq- 
neighborhood of in C™. When H comes from a time-dependent quasi-periodic perturba¬ 
tion of a convex integrable Hamiltonian, the perturbation is in fact independent of J so the 
above requirement is void. Then, if (j^ is satisfied but only for a certain domain B properly 
contained in M™, the statement remains true up to the following modification: the stability 
time T(e) should be replaced by min{To,T(e)} where Tq is the first (possibly infinite) exit 
time of J{t) from the domain Ur^B. This will follow easily from the proof of Theorem [H and 
Theorem [1] is actually a special case as Tq = oo when B = M”*. 
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1.3 Strategy of the proof and plan of the paper 

There are two known methods to prove Nekhoroshev type estimates for small perturbations 
of integrable Hamiltonian systems. The first one is the Nekhoroshev-Poschel’s method, intro¬ 
duced in the seminal work of Nekhoroshev ( |Nek77| . [Nek79j 1 and later improved by Poschel 
l |Pos93j i in the convex case (see also |GCB14j for a further extension of the work of Nekhoro¬ 
shev and Poschel leading to an improved value of the stability exponents in the steep case, 
generalizing the known values of the convex case). The second method is the Lochak method, 
introduced by Lochak l [Loc92| . |LN92| l in the convex case (see also [BN12] for an extension to 
the steep case, though with worse values for the stability exponents). In the convex case, the 
latter method is undoubtedly the simplest and most elegant way to prove stability estimates. 

The Lochak method crucially relies on the existence of periodic orbits for the integrable 
system, that is on the existence of periodic frequencies. Now in the case of a quasi-periodic 
perturbation, the space of frequencies is of the form (a;, a) G where w € M”" is free but 

a G fixed and non-resonant (in particular, a is not periodic; if it were, one would be in 
fact looking at a time-dependent periodic perturbation). The issue is that this space does 
not contain periodic frequencies. As a matter of fact, it is not really necessary to have exact 
periodic frequencies, but only frequencies which can be approximated by periodic ones. Now 
any frequency of the form {uj,a) can be approximated by a periodic frequency, say (a;', a'). 
But then necessarily a' ^ a, and since our Hamiltonian is not convex in the J variables we 
were not able to prove stability close to such periodic frequencies. Therefore the method of 
Lochak does not seem to extend in a easy way to the case of a quasi-periodic perturbation. 

Our strategy is therefore to try to extend the Nekhoroshev-Poschel method, and we will 
succeed in doing so. This method consists of covering the space of frequencies by resonant 
blocks, which are neighborhood of resonances (defined by a certain lattice of integer vectors) 
that are otherwise non-resonant (for integer vectors not in the lattice). Using convexity, 
one can then show that the solutions stay in the same resonant block for a long time (in the 
general steep case, they may leave their resonant block and the proof of the stability gets much 
harder). In the autonomous case where the frequency space is just oj G resonances define 
linear subspaces which are orthogonal to arbitrary submodules of Z"". In the quasi-periodic 
case, our frequency space is still n-dimensional but resonances are associated to submodules 
of Z”"*""*. But not all submodules of are associated to resonances: those that are 

not will be called non admissible. Now resonances associated to admissible submodules do 
not necessarily define linear subspaces but rather affine subspaces in the space of a; G M”. 
Therefore we are facing much more resonances than in the autonomous case, and the geometry 
of these resonances gets more involved. In particular, different admissible submodules might 
lead to different but parallel affine subspaces, and it is at this point that the assumption that 
a is Diophantine (or in fact simply non-resonant) comes into play: it ensures that we can 
control the distance between these parallel affine subspaces. In particular, in the extreme case 
where these parallel affine subspaces are just points (that is, their associated vector space is 
trivial), they can get very close to each other, and this is precisely from this phenomenon 
that our values of the stability exponents come from. This improved geometry of resonances 
will lead to the fact that any frequency (cu, a) is close to some resonant frequency of the form 
(a;',Q:). Using this, and the fact that our integrable Hamiltonian is convex in the I variables 
while linear in the J variables, we will be able to prove stability for the I variables (but we 
will loose control on the evolution of the J variables). 

Now the plan of the paper is as follows. Section ^ deals in details with the geometry of 
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resonances that was alluded above. This section contains the main technical part of the work. 
Section ^ deals with the analysis (construction of a normal form) and the local stability 
results. The analysis, and therefore the stability in the non-resonant case, is completely 
standard and we can simply refer to [Pos93] . The stability in the resonant case uses convexity: 
our integrable Hamiltonian is just “partially” convex so we need to justify that the arguments 
go through, at the expense of loosing control on the J variables. The proof of our main result 
Theorem [1] will be given in Section ^ using the results of Section ^ and Section ^ The 
next sections contain further results that were mentioned in the Introduction: namely, we 
prove better stability results for solutions close to resonances in ^or far way from resonances 
in ^ while in ^we give a more general result assuming a to be only non-resonant. The last 
Section ^consists of concluding remarks. 

2 Geometry of resonances 

The purpose of this section is to study the resonant and non-resonant properties of the 
frequency space 

where a € is a fixed vector, which will be assumed to be Diophantine, and w is a vector 
varying freely in More precisely, our aim is to cover this space by neighborhoods of 
resonances (associated to certain submodules A of on which non-resonant estimates 

can be established (for integer vectors k ^ A). 

2.1 Admissible resonant zones and resonant blocks 

We fix a real parameter K > 1, and in this section, a € will be assumed to be simply 
non-resonant. 

A submodule A of is said to be a AT-submodule if it is generated by elements 

{k,l) G Z"' X Z™ = such that \{k,l)\ < K, and it is said to be maximal if it is 

not properly contained in any other submodule of the same dimension. Given an integer 
1 < d < n -1- m, the set of all maximal K-submodules A of Z"'”*'™' of rank d will be denoted by 
Mfc,d- For A G Mx,d, we define the space of A-resonances by 

i?A = {w G M” I {k, 1) • (cj, a) = 0, V(A:, 1) G A}. (8) 

Quite obviously, since a G M”* is non-resonant, R\ will be non-empty only for certain max¬ 
imal AT-submodules A. Let us consider the subset ^ of MK,d consisting of admissible 
submodules: they are submodules whose intersection with {0} x Z"* C Z” x Z™ = is 

trivial. Equivalently, given any basis {k^, d),..., {k'^, l'^) for A, the vectors k^,... ,k'^ in Z” 
are linearly independent. It is plain to check that if A is not admissible, then R\ is just the 
empty set. Note also that if A is admissible, its rank is at most n. 

Now consider A G ^ where 1 < d < n. If H : ]^n+m ]^n 

is the canonical projection, 

A := n(A) is a submodule of Z"", of rank d, which generates a real subspace (A) of M” of 
dimension d. It is clear that A is a AT-submodule, but it is not necessarily maximal. The space 
of A-resonances defined in (l8|) is non-empty, it is an affine subspace of M”' whose associated 
vector subspace is the vector subspace (A)-*- orthogonal to (A). 

It is not the space of resonances but rather their neighborhoods that will play a role in 
the construction below. To define them, given A G ^ and its associated submodule A, 
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we first define |A| as the co-volume of A viewed as a lattice in (A). It is the volume of the 
fundamental domain spanned by the vectors of any choice of basis for A: letting A be an 
nx d matrix whose columns form a basis for A, then |A| = V detA*A, and this latter quantity 
is easily seen to be independent of the choice of a basis. It is worth recalling, as it will be 
used, that if Sd{A) denotes all square matrices of size d that can be extracted from A, then 
we have the equality (Cauchy-Binet formula) 


|A| = V deiA^A = / ^ (deti?)^. 

Y B&Si{A) 

Then, we introduce n positive real parameters Ai, A 2 , • • •, A„ and, for 1 < d < n, we define 
the associated resonant zone 


where 


ZA:={a;GM’^| ||a;-i?A|| <(5 a}, ^a ^ 

|A| 


leu —iiA||:= inf ||a; —eu'||. 


We then define the resonant zone of multiplicity d, for l<d<n + l, by 


Zd-.= IJ 1 < d < n, Z, 


n+l 


AeM' 


K,d 


The resonant block associated to A G 1 < d < n, are defined by 

Ba ■= Zt^ \ Zd+i 


(9) 


and eventually the resonant block of multiplicity d, for 1 < d < n, is 

Bd := U Ba. 

Setting B{o} = M"" \ Zi, we arrive at the following decomposition 

= B{o} AZi = B{o} U U ^2 = • • • = B{o} U U ... Bn-i U Bn (10) 

since Bn = Zn as Zn+i = 0. 


2.2 Non-resonant domains in freqnency space 

Consider a domain B C M”", a submodule 

AGMf,:= U 

l<d<n 

and a real parameter /3 > 0. Then the domain B is said to be (/3, A')-non resonant modulo A 
if for any {k, I) G such that | (A:, 01 < -A and {k, 1) ^ A and any cu G B, we have 

\{k,l) ■ (w,a)| > 13. 

Our purpose here is to show that the resonant blocks Ba, for A G are (/3A,i^)-non 

resonant modulo A, for a suitable /3 a provided that a is Diophantine and the parameters 
Ai,A 2 ,...,An satisfies certain compatibility conditions. This is the content of the lemma 
below. 
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Lemma 1. Let E > 0 and F > E + 1. Assume that a E M”* satisfies (Dio^^r) o-nd 

IFKXd < Xd+i < l{d + 1 < d < n - 1, 

[An < F-^-i{n + 

Then for any A E M'fi-, the block B\ is {/3 a, K)-non resonant modulo A with 

I3a = EK5a, Ay^{ 0 }, /3|o} = Ai. 


( 11 ) 


( 12 ) 


Proof. Let A E of rank d, with 0 < d < n, and {k,l) ^ A such that \{k,l)\ < K. Let 
A_|_ be the submodule of Z”'"*'™' generated by A and {k,l). Since A is maximal and does not 
contain {k, 1), the rank of A+ is equal to d + 1. 

Let us start with the special case d = 0, that is A = {0}, and fix w E -Bjo}- Either A+ is 
admissible, or not. In the second situation, A: = 0 E Z” while / 7^ 0 E Z”* as {k, 1) ^ A = {0}, 
and we have, using the fact that a satisfies (Dio..y^T-) and m, 

\{k,l) ■ (w,Q!)| = |E q;| > > Ai. (13) 


In the first situation where /c 7 ^ 0 E Z"’, let E Ra+ such that 

llw - w+ll = ||w - Ba+II > dA+. 

Now as {k, 1) ■ {uj+,a) = 0 we have 

{k, 1) ■ (cj, a) = {k, 1) ■ (w, a) — {k, 1) ■ (u;+, a) = k ■ {to — uj+). 

But the vector uj — a;+ belongs to the line orthogonal to (A_|_)-’- = (k)-^ in (A)-*- = M" which 
is nothing but the line generated by k, so we obtain 

\{k,l) • {oj,a)\ = \k ■ (w - w+)| = ||A:||||a; - w+ll > ||A:||(5 a+ = P|||A+|“^Ai > Ai (14) 

where we used the fact |A+| < ||A:||. From (fT^ and (fHll the statement in the case d = 0 
follows. 

Now assume 1 < d < n. It is enough to prove that given any uj E Ra \ Z^+i (where we 
recall that Z^+i = 0 ), we have 

Indeed, for any u E Ba = Za 
llw — d;|| < (5 a and thus 

\{k,l) ■ (d),Q;)| 


\{k,l)-{uj,a)\>FK6A. (15) 

\ Zd+i, by definition there exists uj E Ra \ Zd+i such that 

> \{k,l) ■ {uj,a)\-\k ■ {u - uj)\ 

> \ik,l)-{uj,a)\-\\k\\\\{Cb-uj)\\ 

> \{k,l) ■ {uj,a)\-K\\{uj-uj)\\ 

> FK6a - K6a = {F- 1)K6a > EK6a- 


So it suffices to prove (fT5]) . As before, there are two possibilities for A+: either it is admissible, 
or not. In the second situation (which is obviously the only possibility for d = n), let us choose 
a basis {k^,fi ),..., {k^,l^) for A such that \{kfiP)\ < AT for 1 < j < d. The assumption that 
A+ is not admissible means that A: is a linear combination of A;^,..., k^. For 1 < j < d, let us 
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define E by selecting the first d components of y, and we define k the same way. As 
k},,k'^ are linearly independent (because A is admissible), the determinant A.(k^, ..., k'^) 
of the square matrix of size d whose columns are given by k^,... ,k'^ can be assumed to be 
non-zero without loss of generality. By Cramer’s rule, it then follows that k can be written 
as 

k = vik^ + V2k^ H-+ -+ 

where A = A(fc^,..., fc'’*), Ai = A(fe, ..., k'^), Ad = A{k^, k"^ ... ,k) and for 2 < j < d — 1, 
Aj = A(k ^,..., k^~^, k, ..., k'^). Let us denote 

1^ = 1- - i,dl‘^ e 


The vector is non-zero; if it were, {k, 1) would be a linear combination of {k^,l ^),..., {k'^, k^), 
that is {k, 1) E (A) nZ"'+”*, but since A is maximal it is equal to (A) nZ"'+™' and so this would 
contradict {k, 1) ^ A. Now using the fact that , ) E A for 1 < j < d and uj E R\ we can 

write 

d 

{k, 1) ■ {uj, a) = {k, 1) ■ (w, a) — ^ i'j{k^ ,P) ■ {uj, a) = (/*, a). 

i=i 

The vector AR = Al — Ai/^ — ... — Adl'^ E Z™, and by Hadamard’s formula, we have 

|A| < \\k^\\---\\k^\\ <K^ 
and similarly |Aj| < K'^ for 1 < j < d, so as a consequence 

|A/,| < |A||/| + |Ai||/^| + --- + |Arf||/'^l < (d+l)A:^+^ 


Since a satishes (Dio.y^T-), it follows that 

|A||(fe,Z) • {uj,a)\ = |A||(/„a)| = |(A/„a)| > ^{d + l)-^ 
and as a consequence of Cauchy-Binet formula, |A| < |A| and hence 


Using this last inequality together with m one arrives at 

\{k,l) • (a;,a)| > |A|-iArf+i > |A|-^FiLArf = FK6a. 


(16) 


It remains to treat the case where A_|_ is admissible. Let E i?A+ such that 

llw - w+ll = ||w - d?A+|| > 5 a+- 

Then since {k, 1) ■ (a;+, a) = 0 we have 


{k, 1) ■ {uj, a) = k ■ {uj — w+). 

The vector uj — a;+ belongs to the line orthogonal to (A+)“*“ within (A)-*-, that is it belongs to 
L := (A+)-*"*- n (A)-*- = (A+) n (A)-*-. Let us decompose k = Pk + (Id — P)k where Pk E L 
and (Id — P)k £ P-^. If follows that 

\{k,l) • (w,Q!)| = \k- {uj- w+)| = \Pk ■ {uj - a;+)| = ||PA:||||w - w+|| > ||PA:||(5 a+. 
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But L is also the line orthogonal to (A) within (A+), and therefore |A_|_| < ||PA:|||A| so, 
using m again, 

\ik,l) ■ {oj,a)\ > |A+||A|-i<5a^ = |A|-iAd+i > \A\-^FKXd = FK6 a. (17) 

From (fTHI) and (fT7|) it follows that (fT5]) holds true for 1 < d < n, and this concludes the proof. 

□ 


2.3 Covering by non-resonant domains in action space 

In section ^2.1l we obtained a covering of the frequency space {(w, a) G ~ {a; G M”} by 

resonant blocks associated to admissible submodules, and in Lemma [T] we proved that these 
resonant blocks satisfy some non-resonant properties. 

Recall that we are given an integrable Hamiltonian h which is real-analytic on the domain 
Vr^D, and satisfies (|Mp and (jmj). A subset of D is said to be (/3, A')-non resonant modulo A 
for h if its image by the frequency map V/i is (/3,iL)-non resonant modulo A. It is said to 
be (5-close to A-resonances for h if the Euclidean distance between its image by V/i and the 
space is smaller than 5. 

Pulling back the covering (Hup back to action space using the gradient map V/i and 
using Lemma [U the following proposition will be easily obtained by carefully choosing the 
parameters A^, for I < d < n. 

Proposition 1. For A G of rank d, with 0 < d < n, let us define 

__7_ 

' |A|F"“‘^+l(n -I- l)Ti;^(n+l)T+n-d-|-l ’ 

where we set, by convention, \A\ = 1 if A = {0}. Then there exists a covering of D by subsets 
Da, where A G such that each Da is (/3a, K)-non resonant modulo A and, for A {0}, 
dA-close to A-resonances with 


9MKrA „ mrA 

/3a = «A = —. 

Exactly as in |Pos93j , the introduction of the parameters M and m in the above statement 
is unnecessary (the above proposition does not depend on the assumptions (IMp and (jmj)); 
these parameters are just here for later convenience. 

Proof. Recall that (flOP gives the decomposition 

= R{o} U Hi U ... Bn-l U Bn 


that can also be written as 

U Ba. 

AeM^ 

We now define those resonant blocks by choosing the parameters 


Ad := 


777-7 


F := lOM/m 1 < d < n. 
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With this choice, the inequalities dill) are satisfied (since m < 1), hence Lemma [T] can be 
applied with E := 9M/m < F — 1 = lOM/m — 1 (since M/m > 1). By definition, for each 
non-trivial A of rank d, the block Ba is Ja-cIosc to A-resonances with 


Ad ^_nw_ 

|A| 8F’^-'^+l|A|(n + l)ri;i'{n+l)r+n-d+l 


mrA 


Moreover, using (|12p from Lemma [U is {j5\,K)-non resonant modulo A with 


Pk = EK6a = 


9MKr\ 

8 


while, for A = {0}, is (/3{o}, Ar)-non resonant with 

my ^ 9MK'j 


/3{o} > Ai = 


> 


8F^{n+iyK(^M)r+n - gj^n+l(^^_^iyj^{n+l)T+n+l 

as F > 9M/m. If we define 


9MKr^Qy 

8 


Da:={IgD\ Vh{I) E Ba}, a E M^, 


this defines a covering (up to removing such sets which are empty) of D with all the required 
properties, and this concludes the proof. □ 


3 Normal form and stability estimates 


3.1 Normal form 


Let us come back to our original Hamiltonian (jH|), and recall that the extended integrable 
Hamiltonian is given by 

h{I,J) = h{I) + a-J, {I,J) £ D = D xR^. 

Let us fix A E Quite obviously, if a subset C D is (/3, iL)-non resonant modulo A for 
h, then Z)* := H* x 1^™^ C Zl is (/3, ZL)-non resonant modulo A for h. 

Now a Hamiltonian of the form 




is said to be in A-resonant normal form if 


g{e,ip,I,J)= Y. 

{k,l)eA 


Such Hamiltonians have additional first integrals: indeed, given any vector {uj,a) which 
belongs to the real subspace orthogonal to A, and any solution I{t), J{t)) of the 

system associated to h + g, we have 


(w,a) • 


-(uj,a) ■ d( 0 ^^)(h(I(t),J(t)) + g(0(t),(p(t),I(t),J(t))) 
-(uj,a) ■ d(g^^j(g(0(t), (p(t), I(t), J(t))) 

{k,l)£A 

-z Y {u^,a) ■ 

(A:,Z)gA 

0 . 
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In the special case where A = {0}, it is straightforward to see that g is in fact independent 
of the angles {6, ip), so that the resonant normal form is integrable. 

We can now state the normal form lemma, which states that on a sufficiently small neigh¬ 
borhood of a non-resonant domain modulo A, up to a real-analytic symplectic transformation 
which is close to the identity, the original Hamiltonian can be written as a A-resonant normal 
form h -|- <7 up to an exponentially small remainder. 

Lemma 2. Let H be as in dSl), with h satisfying (IMD and f satisfying m- Consider a 
domain l)* = x C D whieh is {ft, K)-non resonant modulo A for h, and given some 
parameter r > 0 , assume that 


e < 


/3r 


r < 


8/3 


r < ro. 


279A:’ ' “ 9MK' 

and Ksq > 6 . Then there exists a real-analytic symplectic embedding 


(18) 


$ : VfD, X ^ VrD, X VsJ^+^, f := r/2, sq 


so/6, 


such that 

Ho^ = h-\-g-irf* 

where h-\- g is in A-resonant normal form with the estimates 

\g + f*\f,so < 2e, \f*\f,so < (19) 

and 

sup \\Iii^Ml,J,9,ip) -{I,J)\\< ^ ( 20 ) 

where n/,j denotes the projection onto the action space coordinates. 

This statement is a direct consequence of the Normal Form Lemma of |Pos93j (with the 
choice of the constants p = 9/8 and g = 9), to which we refer for a proof. More detailed 
estimates on g and on T are available, but they will not be needed. 

3.2 Non-resonant stability estimates 

In the special case where the domain D* is (/3, K)-non resonant modulo A for h, with A = {0}, 
the normal form obtained in the previous section is, as we already said, integrable up to an 
exponentially small remainder. It is very easy to prove that in this case the action variables 
{I{t), J{t)) remain stable for an exponentially long interval of time, and this does not require 
any convexity assumptions on h. Here’s a precise statement. 

Proposition 2. Let H be as in (0, with h satisfying (IMI) and f satisfying (|^ . Consider a 
domain x C D which is (/3, K)-non resonant modulo A = {0} for h, and given 

some parameter r > 0 , assume that dUl) is satisfied. Then for every solution with initial 
action {Iq, Jq) G I)* we have 

||(/((),J(())-(/„.J„)||<r. |t| < 

This is the content of Proposition 1 (Nonresonant stability estimate) of |Pos93j . to which 
we refer once again for a proof. 
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3.3 Resonant stability estimates 

Next we study the case where the domain D* is (/3, K)-ii.on resonant modulo A for h, with 
A non-trivial. The domain Z)* = D* x will be said to be 5-close to A-resonances if D* is 
5-close to A-resonances, as defined previously. 

Assuming convexity of h, one knows how to bound the variation of the action J{t)) 
using conservation of the energy and convexity arguments, as was first proved in |BG 86 j and 
later in |Loc92j . However, in our situation h is not convex, but it is convex with respect to 
the I variables and linear in the J variables, and we will prove below that this is sufficient to 
bound the variation of the action variables I{t) only, loosing control on the evolution of J{t). 

Proposition 3. Let H he as in (0), with h satisfying (IMp , (jH]) and frmp and f satisfying (jg) • 
Consider a domain = D* x C D which is {P, K)-non resonant modulo A for h, with 
A ^ {0}, but 6-close to A-resonances. Given some parameter r > 0, assume that 


mr'^ mr 8/3 

« = !-' r<r„. 


( 21 ) 


Then for every solution with initial action {Iq,Jq) E D* we have 


P(i) - ^oll < r, 


|i| < 


KsolG 

288Lle 


Proof. First we assume that Ksq > 6 . Since m < M, one easily check that (| 2 ip implies (|18p . 
and therefore Lemma [2] can be applied: there exists a real-analytic symplectic embedding 


^ ^ ^r^soL)^ 


such that 

Ho^ = h + g-\-f^ 


where /i -|- 5 is in A-resonant normal form with the estimates (1191) and 
and ( 1 ^ we get 


||n/,jd>-id|| < 


18K£ ^ 2^e ^ mr 

13 -W?- Wm 


5 


(1^ . Now, from (1^ 


Therefore the inverse image of Z)* x by ‘h is contained in UpD^ x where p := 

5/(4M). Recall that f = r/2 and sq = so/ 6 . We claim that for any initial action (Zq, Jo) E 
UpD^, the solution {I(t), J{t),6(t),(p{t)) of the system associated to the Hamiltonian ZZ o <f> 
satisfies 


||J(t) - Joll <r- P, 


|i| < 


^spr^ Ks/6 

12SJe 


Assuming this claim, for any initial action (Zq, Jq) E D^, the solution (Z(t), J{t), 9{t), ^{t)) of 
the system associated to the Hamiltonian ZZ satisfies 


||J(i) - Joll < ||J(0 - J(i)|| + ||J(i) - ^11 + ||Jo - Joll < p/2 + {r - p) + p/2 < r < r 


for times 

- 2880e 

which is exactly the statement we want to prove. It is therefore sufficient to prove the above 
claim. 
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To simplify notations, let us drop the tildes and simply write {Io,Jo) E C/p-D* and 
J{t),9(t),ip{t)) the associated solution. Let B be the ball of radius r — p around Iq, 
then B x M"* is contained in UfD^. Let Tg be the positive time (possibly infinite) of first exit 
of (/(t), J{t)) from B x R”*: it is then also the time of first exit of I{t) from B. Let also 

r. T:=mm{TM. 

Furthermore, let us write 


Ah := h{I{T), J{T)) - h{h, Jo) E R, 

A(/, J) := (/(T), J(T)) - (Jo, Jo) = (/(T) - Jo, J{T) - Jo) E 
AI := I{T) -IqG R”, 

I{s) := Iq + sAI E R”, 0 < s < 1. 

By definition of h, we have Vh{I,J) = (V/i(/),a) E R”'+™' and 

V^h{I,J) = o) G Mn+mirn, G M„(R). 

Using Taylor’s formula with integral remainder at the point (Jo, Jo) and the special form of 
V^h, we get, letting oj := V/i(/o), 

Ah = (a;,a) • A(/, J) + [ {1 - s)V^h{I{s))AI ■ Aids. 

Jo 

Using the assumption firrii) . that is the convexity of h, we obtain 

|Ah|+ |(a;,a) • A(/, J)| >m/2||A/||2. (22) 

By conservation of energy H and using the first part of (|19p together with the first part 
of (1211) and the dehnition of f, we have 

2 "2 

— mv TUT 

IAh| <2\g + UWs, < 4e < ^ = ^. (23) 

Then, by definition, (Jo, Jo) G UpD^, so there exist (Iq, Jq) G I)* which is p-close to (Iq, Jq), 
but then by assumption Z)* is (5-close to A-resonances, hence there exists w* E R"" such that 
the vector (w*,a) belongs to the real subspace orthogonal to A and is (5-close to (V/i(/q), o;). 
Therefore, since h satisfies ()MD and by definition of /?, we obtain 

||a;-a;*|| = ||Vh(/o) - w,|| = ||Vh(/o) - Vh(/')|| + ||Vh(/') - a;,|| < Mp + 6 = 56/A. (24) 

Writing 

{uj,a) ■ A{I, J) = (cj, a) • A(/, J) — (cj*, q;) • A(/, J)-I-(w*, a) • A(/, J) 

= (cj — cj*, 0) • A{I, J) -I- (w*, q:) • A(/, J) 

= (cj — cj*) • AI + (w*, a) ■ A(/, J) 
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we obtain 


\{LO,a)-A{I,J)\ <\{uj-uj,)-AI\ + \{uj„a)-A{I,J)\. (25) 

Using (I24p . the second part of (I2ip and the definition of f, we can bound the first summand 
by 

7T? ^7T?7^^ T7? 

|(a;-a;.)-A/| < ||a;-u;.||||A/|| < (5V4)||A/|| < _ + -||A/||2 = ^ + _||A/||2. (26) 

For the second summand, using the fact that (a;*,a) belongs to the real subspace orthogonal 
to A, and that h + g' is in A-resonant normal form, we have 

|(a;*,a) • A(/, J)| < f \{ui^,a) ■ 

Jo 

< T\\{ui^,a)\\ sup \\de,^Ml,J,9,ip)\\. 

(I,J,e,ip)eBxR^ xT"+’" 


Now using a Cauchy estimate and the second part of (USD we get 

Pe,^/*(/, J,0,(/p)l| < -^\Mr,s < 

eso eso 

Moreover, asT <T^, and ||(a;*,a)|| < O since h satisfies ([HD) we get 

~2 

|(w*,a)-A(/,J)| (27) 

eso 12e 

Putting together (1^ . ([23D, dSSD, and (f271l we eventually arrive at 


m/3||A/|p < 




+ 


5mf‘^ 

32 


+ 


mr 

12 ^ 


< 


mr^ 


Now p = 6/{AM) = mr/{32M) = mf/{l6M) < r/16 as m/M < 1, thus f — p > 15r/16 and 
in particular < 4(r — p)‘^/3. This, together with the last inequality, implies that 


m/3||A/|p < m/3{r — p)^ 


and therefore ||A/|| < f — p. This eventually proves that T = T^, that is 


ll-f(t) - -foil <r- p, 


0 < t < T* = 


mspr^ Ks/o 

l2Q.e 


The same argument yields the same result for negative times —T* < t < 0, and hence 


\\I{t) - -foil <r- p, 


|f| < 


mspr^ Ks/6 

l2Q.e 


which was the claim that needed to be proved. This ends the proof under the assumption 
Ksp > 6. But if Also < 6, the exact same argument applies to the original Hamiltonian 77, 
by setting g = 0 and f* = f, and this concludes the proof. □ 
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4 Proof of the main result 


This section is devoted to the proof of our main result, Theorem [H which will be easily 
obtained using Proposition [H Proposition [2] and Proposition [3l 


Proof of Theorem [7J Recall that we are considering H as in (iHli , with h satisfying (IMh , (jH]) 

. Recall also 

(n + l)r + 1 


and fiTrlj) , / satisfying (j^ , and a satisfying (Dio^,,-). Recall also that we have defined the 
positive constants a and b by 


1 


a = 


b = 


2 (n + l)(r + l)’ “ 2 (n + l)(r + l)' 

We can apply Proposition [T] to obtain a covering of D = D x by subsets D\ = D\ x M™', 
where A G such that each D\ is (/3a, R)-non resonant modulo A for h and, for A 7^ {0}, 
(^A-close to A-resonances for h with 

9MKr\ 


/3a = 


where 


TA : = 


8 

7 


(3a = 


mri^ 


d = rank A. 


(28) 


(29) 


_j_ ly J^{n+l)T+n-d+l ’ 

Then, given any A G Mf., we can apply either Proposition [2] (for A = {0}) of Proposition [3] 
(for A / { 0 }), with r = ta, /3 = /3a and 6 = 5\ provided that (fT8|) and (f2T|) are satisfied. The 
inequalities (fT8]l are easily seen to be implied by (1211) as we already pointed out, hence we 
only need to verify m, and in view of our definitions of r\, /3a and 5a, the latter reduces to 


e < 


mr\ 


ta < ro- 


(30) 


Given any A G we have 


7 


<rA< 


7 


_pn+i(„ I)r|;^(n+1)(T+1) - “ F(n + 1 )'^A:("+i)^+i ’ 

as F > 1 and K > 1. Hence (15UD is satished, for any A G if 


£ < 


my 


2l0_p2(n+l)(^ lyr J^2(n+l){T+l) ’ 


K > 


7 


Let us define 


and 


eo ;= 


my 


2l0^2(n+l)(„ 1)2 t’ 


e* So 


F{n + l)^ro 

F{n + l)^ro^ ^ 
7 


(31) 


R:= 


then ([3T]) eventually reduces to 

e < min{eo,e*}- (32) 

Under this choice of K and this smallness assumption on e, Proposition [2] and Proposition [3] 
apply and in the resonant case, we obtain 

b 


II^W-^o||< 


7 


< 


7 


F(n + l)^Ar(’^+iR+i - F(n + 1) 


£0 
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up to times 


1^1 < T\ exp 
while in the non-resonant case, we obtain 


£0 /£o\“ 
6 V e 


, Ta: = 


msorj^ 
2880e ’ 


||/(t)-/o||<||(/(i),J(t))-(/o,Jo)||< 


< 


F(n-I-F(n-|-1)"^ \eo 


up to times 


So /eo\“ 


t| <T{o}exp(-^(^-j ], r^o} := 


SQ^’jO} 

5e 


To obtain a uniform time estimates, observe that for any A G we have 


288 o “ n 


Letting 


we have just proved that 


R* = 


7 


F{n + 1)' 


rj. _ 3S0 

* O’ 


eo 


||L(t) - loll < — , |t|<T*exp — 


So (eo 


6 V e 


provided (|32p is satisfied. Recalling that F = lOM/m, this was exactly the statement to be 
proved. □ 


5 Improved stability close to resonances 

For solutions starting close to resonances, we can obtain a better result. Consider a fixed 
submodule L of of rank d, which is assumed to be admissible and maximal, and let 

Kl > 1 such that L is a iF^-submodule. Recall that L denotes the projection of L onto EF, 
and |L| denotes the co-volume of L. To such a L the space of L-resonances is 

Rl = {w G M” I {k,l) • (w,a) = 0, 'i{k,l) G L} 


and let 

Sl := {(/, J) gD \ V/i(/) £Rl} = {I €D \ Vh{I) G Rl} x 

the resonant domain in action space. Solutions with initial action close to satisfy better 
stability estimates, as stated in the theorem below. 


Theorem 2. Let H be as in dSl), with h satisfying (|MD . ([11]) and ([Trljl . and f satisfying (Q . 
Assume also that a satisfies (Dio.y^T-); o,nd let us define 

1 


a{d) = 


2 (n -|- l)r -|- n -|- 1 — d’ 


m = 


{n -k 1 )t -k 1 


2 (n -|- l)r -|- n -|- 1 — d’ 



n-\-l—d 


7 

(n + 1)^’ 
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and, given a submodule L as above, let 


£o{L) 


1 m 7 ^ f m N^2(n+l-d) fr\ f ^0 ^ 

^2i0(n + l)2- IiomJ ’ ^ J 


e**(L) 

//e < mm{eo{L),e,t{L),£^:,f{L)}, for any solution {I{t), J{t),9{t),ip{t)) of the system asso¬ 
ciated to H with initial condition {Iq, Jo,6o,ipo) G UpSi x with p = 4:M~^y/me, we 

have 

||/W-/oil <«.(<!) (^) , |t|<r.exp(^|(2!S) j 

with T* as in Theorem [7J 

Observe that in the special case where L = {0}, Sl = Dx M™, \L\ = Kl = 1 and therefore 
the above statement exactly reduces to Theorem [TJ 

Proof. Consider the subset 


M^(L) := {A € I A D L} 


and assume, for any A G Mf^{L), 

rA<ro, K>Ki^ (33) 

where the last inequality is to ensure that Mf^{L) is actually non empty. 

The resonant blocks i?A, for A G Mf^{L), cover the resonant zone Zl and hence their 
pullbacks cover UpSi, with 


, 6l mri ^ I - ! - 

p = — =- > - Pme = 4m Jme = 


Moreover, for 


P- 


_7_^_7_ 

\l\pn-d+l(^^ -^yj,.(n+l)T+n-d+l - Ta _ ^ • 

Using the above inequalities and proceeding exactly as in the proof of Theorem [U we can 
define 

A':= (""Wy"' 

and verify that ([551) is implied by 


£<£o{L), £<£^{L), £< (L). 


The stability estimates apply uniformly to all blocks Ba, for A G and one easily 

check that the statement follows with the given constants. □ 
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6 Improved stability far away from resonances 

We now investigate solutions which start far away from resonances. In this special case, we 
will actually be able to control the variation not only of the I but also of the J variables. 

Results of this section do not depend on h being convex or its gradient being bounded, 
that is l[r7i]l and m are unnecessary. It will be sufficient to assume the existence of m > 0 
such that for any Lebesgue measurable subset H G M"', we have the measure estimate 

Leb(V/i~^(n) r\ D) < ?fi~^Leb(n) {fh) 


where Leb denotes the Lebesgue measure on M”. Certainly, the convexity assumption (fm^ 
implies (I ml) with fh = m, but the latter is more general: in particular, it holds true if h is 
Kolmogorov non-degenerate, that is if the determinant of the Hessian Vh?‘{I) is uniformly 
bounded away from zero for any I ^ D. 

Now consider A € admissible maximal iL-submodule of of rank 1. Then 

A contains a unique vector {k,l) G Z"’ \ {0} x Z™ such that |(A:,/)| < K and such that 
its components are relatively primes. The submodule A of Z"" is then generated by k, and 
|A| = ||A:||. In the sequel, we shall write A = A{k,l). 

With these notations, we recall that the completely non-resonant block H|o} C M”, intro¬ 
duced in Section ^ can be defined by 


-^{0} s w G 


— RA{k,l)\\ > 


Ai 

ll^ll 


Ai — 


7 


8F^{n + ’ 


F = IQM/m. 


Therefore, from Proposition [3l this set is (/3{o}, Ar)-non resonant with 


/3{0} 


dMKrfQ'x 


_7_ 

pn+i(^n -I- l)'rAr(’^+i)C+b ’ 


If we further define 


-^{0} ■— {I ^ D \ Vh{I) G -^{0} T>{o} X M"*, 


we arrive at the following statement. 

Theorem 3. Let H he as in dSl), with h satisfying (IMI) and (Iml) . and f satisfying (jg). 
Assume also that a satisfies (Dio..y^T-); and let us define 

n +1 


R* ;= 


( m \ 
VIOM/ 


7 


(n + 1)^’ 


£0 


MY 


2i0(n-h 1)2^ VIOM 


/ m '^2(n-|-l) 

V; 


) 


. ^0 
= £0 I ^ 

-it* 


If e < min{eo,£*}, for any solution (/(t), J(t), 0(t), of the system associated to H with 
initial condition {Iq, Jo,do,I>o) € x we have 


\\{I{t),J{t))-{Io,Jo)\\<R* 



|i| < 


SqRi^ 

5eo 



Moreover, the complement of iA{o} in has a measure of order , with a and b as in 
Theorem [7J 
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Observe that the measure estimate on the complement of D{o} is better than the one 
obtained in |Pos93j for r = 0, as in the latter reference it is only of order one. Observe also 
that it is only for this measure estimate that ([rn|l is needed; the first part of the statement 
holds true without this assumption. 

Proof. Proposition [2] can be applied with r = r{o}, provided that 


- 2’^9K 


^io} 

210 ^ 


?’{o} < ’’o, 


(34) 


holds true. Choosing 


we have 


K := 



r{o} = R* 



1 

2 


and the inequalities ([54|1 are satisfied if e < min{eo, where eo; and R* are the constants 
given in the statement. It then follows from Proposition [2] that 


\\m,J{t)) - (Jo, Jo)|| < r^o}, |t| < 

that is 





t 


SqR,:^ 

5eo 




which proves the first part of the statement. 

Concerning the second part of the statement, we follow [Pos93j . First, in view of (|ml) . it 
suffices to show that the complement of I?{o} in has a relative Lebesgue measure of order 
£^. But by construction, the latter set is Zi, the resonant zone of multiplicity 1 defined in 
Section ^ which has a relative measure of order 


E 


A(fc,0eM^ 


Ai 


PI I _fi:(n+l)r+n 


E 


0<\k\<K 


PI I /i^(n+l)T+l 


~ e 


a((n+l)T+l) _ 


□ 


Next we look at a different, and in some sense more natural, non-resonant set. Choose 
0 < 7 ' < 7 and t' such that t' > n + m — 1 and t' > r. We define 

By,r' ■■= {wEM’^I \ik,l)-iuj,a)\ > jf\k\ + ^ 


and 


D 


r,T' 


:= {I e D \ Vh{I) E Dyy := 


r' X 


Those sets are clearly (/?', iir)-non resonant, with /3' := ^'K , and we obtain the following 

result. 
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Theorem 4. Let H he as in dHl), with h satisfying (IMIl and (Irnl . and f satisfying (jg). 
Assume also that a satisfies (Dio^^T-j, and let us define 


a := 


1 ' 


,/2 


2(r' + l)’ *■ 9M’ 768M’ 


e* — ^0 




If e < min{eQ, e'^}, for any solution {I{t), J{t),9{t),ip{t)) of the system associated to H with 
initial condition {Iq, jQ,6o,ipQ) G Dyy x we have 


mt),j{t))-{io,Jo)\\<K 


1 

e N 2 


< 


SqK 

5e'n 


0^" 


exp 


/ / \ <1 
£ 0/^0 

6 


Moreover, the complement of Dyy in M” has a measure of order y'. 

As we already said, when m > 2, the set of vectors a G M”* that satisfy (Dio^^,-) with 
T = m — 1 has zero measure, and when r > m — 1, we can choose t' = n + t and thus 

1 


a = 


2(n + r + 1) 


in the statement above. 
Proof. Let us define 


r : = 


8/3' 


87' 


9MK DMiL^'+i ■ 
Proposition [2] can be applied with r = r', provided 

/3'r' Mr'‘^ 


£ < 


2'^9K 210 


r’ < ro 


(35) 


holds true. Choosing 


we have 


K : = 


f \ a 


r' = E!A^ 


1 

e \ 2 


and the inequalities (l3^ are satisfied if e < min{eQ,e'^}, where a', Eq, and ii( are the 
constants given in the statement. Then, exactly as before. Proposition [5] yields 


||(I(t),J(t))-(Io,Jo)||<K -7 


1 

£ \ 2 


|i| < 


sqK (" 

5e' 


exp 


which gives the first part of the statement. 

Concerning the second part of the statement, as before it is enough to prove that the 
complement of By y has a relative measure of order 7'. The complement of By is 

{u; G M- I 3{k,l) G Z-+-, \{k,l) • (a;,a)| < y'{\k\ + |/|)-"'} , 

but since a satisfies (Dio^y^r), and since 7^ < 7 and r' > r, this set is also equal to 

IwGM''! 3{k,l) gZ’"\{0} xZ”', \{k,l)-{uj,a)\ < 7'(|A:| + |, 




But now the above set is known to have a relative measure of order 7': this is exactly the 
content of Lemma 2.12 in [.Tor91| . □ 
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7 A more general stability result 


Let us finally give a more general result, where the Diophantine assumption (Dio^^r) is re¬ 
moved. Assuming a E to be simply non-resonant, we can define a function 'L = 'Lq 
by 

^(iL) = max{|fc-ar^ I fcEZ”^, 0 < |fc| < a:}, K > 1. (36) 

Then we define A = Aq, by 


A(x) = supjiL > 1 I K'1/(K) < x}, X > 'L(l) = \a 


oo 5 II oo • — max I rr d • 

l<j<m 


(37) 


If a satisfies (Dio.y^T -)5 then the functions T and A defined above satisfy 


'L(A') < 7 A(x) > (7x)-r+i. 


(38) 


The only place where (Dio^y^r) was used was in Lemma [TJ But using the function dt instead, 
the exact same proof yields the following more general lemma. 


Lemma 3. Let E > and F > E + 1. Assume that a E 


is non-resonant and 


I FKXd < Xd+i < ^{{d + l)iL"'+i)-i, 1 < d < n - 1, 

[An < F-iT((n + 1 )K^+^)-^. 

Then for any A E the block B\ is {13 a, K)-non resonant modulo A with 


f3A = EKdA, A^{0}, /3|o} = Ai. 

Using this Lemma instead of Lemma [H we arrive at the following proposition which 
generalizes Proposition [TJ 

Proposition 4. For A E of rank d, with 0 < d < n, let us define 

1 

\A\F^-<i+'^^{{n + 


where we set, by convention, \A\ = 1 if A = {0}. Then there exists a covering of D by subsets 
Da, where A E such that each Da is {Pa, K)-non resonant modulo A and, for A ^ {0}, 
dA-close to A-resonances with 


Pa 


9MKrA 

8 


-5a 


mrA 

~ 8 ~ 


Then, using Proposition 0] instead of Proposition [TJ together with Proposition [2J and 
Proposition [3j and proceeding exactly as in the Proof of Theorem [TJ we obtain the following 
statement. 


Theorem 5. Let H be as in with h satisfying (IM[) . (IHj) and (irTTl) . and f satisfying (jg). 
Assume also that a is non-resonant, and that 

< (n + l)^m|a|^ 

^ - 2^0i?2(n+l) 
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so that we can define 


A, 


/ n + 1 fm\ 
■= ^ ( 25F«+i' ~ 


e 


Ke ■■ = 


A, 


n + 1 


1 

71 + 1 


: = 


FKs^{{n + 1)A:< 




r/ien, if K^>1 and Re < ro, for any solution {I{t),J{t),6{t),if{t)) of the system associated 
to H with initial condition {Iq, Jo,6o,^po) E D x we have 


\\I(t) - IqW < Re, |t|<r*exp 


( spKe 

V 6 


When a is Diophantine, then (l38|) holds true and the above statement exactly reduces to 
Theorem [TJ There also analogues of Theorem [2] and Theorem [3] in the above context that can 
be stated and proved the same way; details are left to the reader. 


8 Concluding remarks 


We conclude this paper with several remarks. We first discuss the Chirikov-Lochak conjecture, 
and then the possibility of extending our results to less regular Hamiltonians, more general 
classes of integrable Hamiltonians and finally more general time-dependence. 

First let us recall that the conjecture of Chirikov-Lochak predicts that 


a = b = 


1 

2(n -I- 1 -I- r) 


(39) 


while we proved 


2(n + l)(l + T)’ ^ 


(n -h l)r 1 
2(n + l)(l+T)’ 


For r = 0, the exponents are the same but not for r > 1. However there is certainly no 
contradiction here; our result yields a better confinement but on a worse time-scale. Chirikov 
initial conjecture is ultimately based on the ansatz 


\{k,l) ■ {ui,a)\ ~ 


1 

^n+m—1 ’ 


|A:| + |/| <iL, 


(40) 


for an arbitrary fixed vector a. E M™, leading to the exponents 


a = b 


1 

2(n -|- m) 


(41) 


We refer to |CV96] . Equation (1.8) for instance. When m = 0 or m = 1, that is when the per¬ 
turbation is autonomous or time-periodic, vectors w E M” satisfying the above condition (140 p 
do exist: they are called badly approximable, they are the “best” non-resonant vectors and 
they form a dense set but of zero Lebesgue measure. In those two cases m = 0 and m = 1, the 
exponents (i4T]l turned out to be correct, but proofs do not use existence (nor density) of those 
badly approximable vectors: on the contrary, in the work of Lochak it is the “worst” resonant 
vectors, namely periodic vectors, that play the major role. When m > 2, as observed by 
Lochak in |LMS03j . the ansatz (1401) and hence the exponents (|41ll are clearly too optimistic: 
one should require a to be Diophantine with some exponent r > m — 1, and replace iOl) by 


\{k,l) ■ (a;,a)| ~ 


1 

J^n+T ’ 


|A:| + |/| < AT, 


(42) 
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so that the exponents of (HT|) become those of ([39]) . Vectors w E M” satisfying ()i2]l do exist, 
they form a dense set which is furthermore of full Lebesgue measure when r > m — 1 (when 
r = m — 1 they have zero Lebesgue measure). In our opinion, it is a very interesting open 
problem to derive first these exponents when m = 0 or m = 1 using badly approximable 
vectors, to see if and how one can try to use vectors satisfying (H2I1 to possibly reach the 
values (f3^ . 

Then, let us comment on the regularity assumption in the results we proved. The assump¬ 
tions that h and / are real-analytic were only used in the normal form Lemma[2l which is taken 
from [Pos93j . So in order to have results for Hamiltonians which are not real-analytic, for 
instance Hamiltonians which are only Gevrey regular or finitely differentiable, one just needs 
a version of Lemma [2] in those settings, and this appears to be only a problem of technical 
nature. One can find in |Boul3bj and [Boul3cj normal form results for non-analytic Hamilto¬ 
nians in the spirit of Lemma [21 even though those statements do not recover Lemma |21 Let us 
also recall that in the autonomous case (or time-periodic case when h is convex), Nekhoroshev 
type estimates are known for Gevrey or finitely differentiable Hamiltonians f [MS02j . [BoulOj . 
|Bonll| l but all those proofs are based on the Lochak method. 

Next let us discuss the more interesting question of whether our results extend to steep 
integrable Hamiltonians, the original class of integrable Hamiltonians considered by Nekhoro¬ 
shev. One should first recall that in the general steep case, unlike what happens in the 
convex case, motions near resonances are not necessarily confined: if they are not, the steep¬ 
ness property ensures that they evolve towards a less resonant domain and, eventually, end 
up in a non-resonant domain on which an integrable normal form (up to a small remain¬ 
der) can be constructed, leading to the stability of the action variables. The issue is that 
for a quasi-periodic perturbation, the non-resonant domain might be empty so the extension 
to steep integrable Hamiltonians does not seem straightforward. Indeed, consider in fact a 
time-periodic perturbation of the simplest class of steep Hamiltonians, namely quasi-convex 
Hamiltonians. Furthermore let’s look at the simplest quasi-convex Hamiltonian, which is 
certainly 

/l(/i, ...,/„) = -{ll H-h In_i) + In- 

After a periodic perturbation the extended integrable Hamiltonian then reads 

h{Il, . . . ,In,J) = 2 + ■ ■ ■ + ^n-l) + In + J 

and frequencies associated to h are all resonant (as the last two components of any frequency 
vectors are both equal to 1). Thus for such a Hamiltonian, integrable normal forms cannot be 
obtained (and even “almost” integrable normal forms, in which the dependence on the angle 
associated to J is allowed, as one can easily check) and to prove stability, if any, one should 
follow a different path than the one of Nekhoroshev. 

To conclude, let us mention that it seems quite unlikely to have a non-trivial stability 
result for an arbitrary time-dependent perturbation, unless the time depends on the small 
parameter, in which case the conjugated action variable can be considered as degenerate (see 
[GZ92] or |Boul3aj i. But it may be possible to extend our results for a class of perturbation 
whose Fourier transform (with respect to time) has suitable localization properties. 
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